Evolution of superconducting order in Pr(Osi_j;Rua;)4Sbi2 
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We report measurements of the magnetic penetration depth A in single crystals of 
Pr(Osi_2,Rua;)4Sbi2 down to 0.1 K. Both A and superfluid density ps exhibit an exponential behav- 
ior for the a;>0.4 samples, going from weak (a;=0.4,0.6), to moderate, coupling (a;=0.8). For the 
a;<0.2 samples, both A and ps vary as T^ at low temperatures, but ps is s- wave-like at intermediate 
to high temperatures. Our data are consistent with a three-phase scenario, where a fully-gapped 
phase at Td undergoes two transitions: first to an unconventional phase at rc2jS^ci, then to a nodal 
low-T phase at Tc3<Tc2, for small values of x. 



The recent discovery [ll |^ of the heavy Fermion 
(HF) skutterudite superconductor (SC) PrOs4Sbi2 has 
attracted much interest due to its differences with the 
other HFSC. Early work suggested that the ninefold de- 
generate J = 4 Hund's rule multiplet of Pr is split 
by the cubic crystal electric field, such that its ground 
state is a nonmagnetic doublet, separated from the 
first excited state Fs by ^ 10 K. Hence its HF behav- 
ior, and consequently the origin of its superconductiv- 
ity, might be attributed to the interaction between the 
electric quadrupolar moments of Pr'^+ and the conduc- 
tion electrons [l| . More recent results appear to rule this 
mechanism out, giving strong evidence for a singlet Fi 
ground state with a F5 triplet state at a slightly higher 
energy [3, 3 . In this scheme, aspherical Coulomb scat- 
tering 01 and spin-fluctuation scattering j5| have been 
proposed as mechanisms leading to superconductivity 

Surprisingly, replacement of Os by Ru, i.e. in 
PrRu4Sbi2, yields a superconductor with K 
and significantly different properties. The effective 
mass of the heavy electrons calculated from de Haas- van 
Alphen (dHvA) and specific-heat measurements P, 0] 
show that, while PrOs4Sbi2 is clearly a HF material, 
PrRu4Sbi2 is at most, a marginal HF. Various experi- 
mental results suggest that these two materials have dif- 
ferent order-parameter symmetry. Firstly, there is no 
Hebel-Slichterpeak in the nuclear quadrupole resonance 
(NQR) data [S| for PrOs4Sbi2, while a distinct coher- 
ence peak was seen in the Sb-NQR 1/Ti data for 
PrRu4Sbi2. Secondly, the low-temperature power-law 
behavior seen in specific heat Q and penetration depth 

S, and the angular variation of thermal conductivity 
, suggest the presence of nodes in the order param- 
eter of PrOs4Sbi2. Specifically, Refs. and 111 reveal 
the presence of point nodes on the Fermi surface (FS). 
For PrRu4Sbi2, however, exponential low-temperature 
behavior was seen in 1/Ti |^ and penetration depth 
data. The latter data were fit with an isotropic zero- 
temperature gap of magnitude A(0)=1.9A:_b?c, showing 
that PrRu4Sbi2 is a moderate-coupling superconduc- 



tor. Thirdly, muon spin rotation (/iSR) experiments 
on PrOs4Sbi2 reveal the spontaneous appearance of 
static internal magnetic fields below Tc, providing ev- 
idence that the superconducting state is a time-reversal- 
symmetry-breaking (TRSB) state J^] . Such experiments 
have not been performed on PrRu4Sbi2. 

It is puzzling that the substitution of Ru for Os (same 
column in the periodic table) causes PrRu4Sbi2 to dif- 
fer in so many respects from PrOs4Sbi2, particularly if 
symmetry of the superconducting gap varies as we go 
from Os to Ru. Recently, Frederick et al. performed x- 
ray powder diffraction, magnetic susceptibility and elec- 
trical resistivity measurements 0| on single crystals of 
Pr(Osi_a;Rua;)4Sbi2. They found a smooth evolution of 
the lattice constant and Tc with x, albeit with a deep 
minimum (0.75 K) in Tc at a;=0.6, and an increased 
splitting between the ground and excited states of the 
Pr ion. On the other hand, one still has to contend 
with measurements |0, 0, 0, 0| that indicate point- 
node gap structure. TRSB and a double superconducting 
transition Tc2<Tc UJ in PrOs4Sbi2, none of which are 
seen for a:>G. We report here a complementary study of 
Pr(Osi_a;Ru;r)4Sbi2 using the penetration depth. 

A recent paper observed an unexpected enhance- 
ment of the lower critical field Hci{T) and the critical 
current Ic{T) deep in the superconducting state below 
T«0.6 K (r/Tc«0.3) in PrOs4Sbi2. They speculate that 
this reflects a transition into another superconducting 
phase that occurs below TcS^O.G K, and may explain 
anomalies in other measurements, such as the levelling 
off of Sb-NQR l/Ti below 0.6 K 0, the smafl downturn 
of penetration depth below 0.62 K and its deviation from 
point-node-T^-behavior above ~0.6 K [loj . 

In this Letter, we present high-precision measurements 
of the penetration depth A(T) of Pr(Osi_2;Rua;)4Sbi2 
(x=0. 1,0. 2,0.4,0.6,0.8) at temperatures down to -0.1 K 
using the same experimental conditions as for PrOs4Sbi2 
and PrRu4Sbi2 10, 12] . For the a::>0.4 samples, both 
A(T) and superfluid density Ps{T) exhibit exponential 
behavior at low temperatures, supporting the presence of 
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an isotropic superconducting gap on the FS. The Ps{T) 
data agree with the theoretical curve over the entire tem- 
perature range. The values of A(0) used in the fits sug- 
gest an increase in coupling strength from weak-coupling 
(a;=0.4,0.6) to moderate coupling (x=0.8). On the other 
hand, the x<0.2 samples exhibit a low-T power law, im- 
plying the existence of low-lying excitations. However, 
the ps data fit a fully-gapped theoretical curve from in- 
termediate temperatures up to Tc, but not curves based 
on a superconducting gap with line or point nodes. This 
is consistent with the scenario depicted by Cichorek et 
al. 13 1 where for the a;<0.2-samples, the fully-gapped 
high-T phase undergoes a transition into a nodal low- 
T phase below Tc3{x). As x increases, the low-T phase 
is suppressed {Tcs decreases) such that for the x > 0.4- 
samples, Tc3 falls below the base temperature of our ex- 
periment, and we are left with a fully-gapped phase over 
our entire experimental temperature range. Taken to- 
gether with other data, we suggest that, in addition to 
the two phases at Td and Tc2 , there is a third supercon- 
ducting phase at Tc3 that exhibits point nodes. 

The single crystal samples were grown by Sb self-flux 
method The observation of dHvA effect both in 

PrOs4Sbi2 and PrRu4Sbi2 could be an indirect evidence 
of high quality of these samples grown in the same man- 
ner. Measurements were performed utilizing a 21-MHz 
tunnel diode oscillator j^l with a noise level of 2 parts in 
10^ and low drift. The magnitude of the ac field is esti- 
mated to be less than 40 mOe. The sample was mounted, 
using a small amount of GE varnish, on a single crystal 
sapphire rod. The other end of the rod is thermally con- 
nected to the mixing chamber of an Oxford Kelvinox 25 
dilution refrigerator. The sample temperature is moni- 
tored using a calibrated Ru02 resistor at low tempera- 
tures (Tbase~l-3 K) and a calibrated Cernox thermome- 
ter at higher temperatures (1.2 K-1.8 K). 

The deviation AA(r)=A(T)-A(0.1 K) is proportional 
to the change in resonant frequency Af{T) of the os- 
cillator, with the proportionality factor G dependent on 
sample and coil geometries. We determine G for a pure 
Al single crystal by fitting the Al data to extreme non- 
local expressions and then adjust for relative sample di- 
mensions [l^ . Testing this approach on a single crystal 
of Pb, we found good agreement with conventional BCS 
expressions. The value of G obtained this way has an 
uncertainty of ±10% because our samples have a rectan- 
gular, rather than square, basal area !l!||. 

We first discuss the x > 0.4 samples. Figure Q (O) 
shows AA(r) for the three samples (a;=0.4,0.6,0.8) as 
a function of temperature in the low-temperature re- 
gion. The insets show AA(r) for the entire temperature 
range. The onset of the superconducting transitions T* 
are 0.81 K (a;=0.6) and 0.88 K(x=0.8). These values 
are consistent with those of Ref. |lj. We could not obtain 
T* for the a;=0.4 sample as the ac losses were so large 
that oscillation was lost before Tc was reached; its large 
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FIG. 1: (O) Low-temperature dependence of AA(T) for (a) 
x^OA, (b) x=0.6, and (c) a;=0.8. Lines: fits to BCS low-T 
expression from Tbase to 0.4Tc. The parameters of the fits 
are described in the text. Insets show AA(T) over the full 
temperature range. 

transition width is also consistent with the ac suscepti- 
bility data of Frederick et al. |Q| , though the origin is 
unknown. The values of Tc, determined from the point 
where the experimental superfluid density almost van- 
ishes and fit the theoretical curves (described later), are 
0.8 K (x=0.4), 0.76 K (x=0.6) and 0.86 K (x=0.8). 

For all three samples the data points flatten out be- 
low 0.3Tc, implying activated behavior in this tem- 
perature range. We fit these data to the BCS low- 
temperature expression in the clean and local limit, 
from Tbase (~0.1 K) to 0.4Tc, using the expression 
AA(T)oc^7rA(0)/2fcBTexp(-A(0)/feBT), with the pro- 
portionality constant and A(0) as parameters. The 
best fits (solid lines) are obtained when A(0)/fcsTc=:1.64 
(x=0.4), 1.53 (x=0.6) and 1.95 (a;=0.8). This implies 
that the a;=0.4 and 0.6 samples are weak-coupling, while 
the 2;=0.8 sample is a moderate-coupling, supercon- 
ductor. The x—0.8 result is consistent with that for 
PrRu4Sbi2 (x=l). 
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TABLE I: Parameters used to calculate curves in Figs. 2 and 
3. Values for x—0 and x=l are included for comparison. 
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FIG. 2: (O) Superfluid density ps{T) = [A^(0)/A^(r)] calcu- 
lated from AA(T) data in Fig. [T] for (a) a;=0.4, (b) a;=0.6, 
and (c) 2;=0.8. Lines; Theoretical Ps{T) with parameters 
A(0)/fc_BTc and AC/'yTc mentioned in the text. 

To extract the superfluid density ps from our data, we 
need to know A(0). Absent published data on A(0), we as- 
sume that it lies in the vicinity of 344 nm (for PrOs4Sbi2) 
113 and 290 nm (for PrRu4Sbi2) 0. We compute 
Ps for an isotropic s-wave superconductor in the clean 
and local limits using ps = 1 + 2 §e^^' where / = 
[exp{E/kBT) + l]^^ is the Fermi function, and E = [e^ + 
A(T)^]^/^ is the quasiparticle energy. The temperature- 
dependence of A(T) can be obtained by using (iH 
A{T)^dsckBT,teiuh{{TT/Ssc)Vi2/m^C)/C][{T,/T)~ 1]}, 
where Ssc=A{0)/kBTc is the only variable parameter. 
The specific heat jump AC/C can be obtained from 
A{0) /ksTc using strong-couphng equations [Ull^l . 

Fig-Elshows the experimental(O) and calculated (solid 
line) values of ps as a function of temperature for the 
X > 0.4 samples. The theoretical curves fit the data very 
well using the parameters shown in TableQ] Fitted values 
for A(0) are reasonable, considering the uncertainty in 
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FIG. 3: (O) Low-temperature AA(r) for (a) a;=0.1 and (b) 
a;=0.2. Lines: fits to AA(T)=A-HBr" from 0.1 K to 0.53 K. 
Insets show AA(r) over the full temperature range. (Q) Su- 
perfluid density Ps{T) calculated from AA(r) data for (c) 
a;=0.1 and (d) 2:=0.2. Lines: Theoretical ps(T) with weak- 
coupling parameters. Note the deviation of data from the 
theoretical curve at low temperatures is more pronounced for 
a;=0.1 than for a:=0.2. 

obtaining the calibration factor G. 

We now turn to the x<0.2-samples. Figs. and|3|3 
show AA(r) in the low-temperature region. The insets 
show AA(r) for the entire temperature range. T* is mea- 
sured to be 1.76 K (a;=0.1) and 1.77 K (2;=0.2), while 
Tc is 1.4 K (2;=0.1) and 1.2 K (x=0.2). A fit of the 
low-temperature data (up to 0.53 Kw0.3T*) to a vari- 
able power law A\{T)^A+BT"- yields 7i=2.5 (x^O.l) 
and 3.3 (x=0.2), indicative of low-lying excitations and 
incompatible with an isotropic gap. 

Figs. ISt and show the experimental (Q) values of 
Ps{T). The solid lines represents the theoretical curve 
based on an isotropic weak-coupling gap as in Table ^ 
Note that the data do not agree with the theoretical curve 
at low temperatures, but agree from intermediate tem- 
peratures up to near Tc- The deviation of data from the 
theoretical curve at low temperatures is more pronounced 
for a;=0.1 than for a;=0.2. This is consistent with the 
scenario depicted by Cichorek et al. 0| , where for these 
low- a; samples, the fully-gapped high-T phase undergoes 
a transition into a nodal low-T phase below Tea (x) . Our 
data also agree with the theory of Hotta , which pre- 
dicts that when the Fi-T .-j sp acing increases (observed as 
X is increased from to 1 |l4|. and for x=l 0), supercon- 
ductivity changes from unconventional to conventional. 
We assume that this nodal phase is a point-node one, 
consistent with Refs. liolllll and so AAocT^ in this phase. 
Consequently, we plot AA(T) vs T^, shown in Fig.^ and 
^p. Tc3{x) is determined from the temperature where 
the data deviate from linearity, from which we obtain 
rc3(a;=0.1)Ri0.32±0.02 K and rc3(x=0.2)=0.15±0.02 K. 
Together with Te3(a;=0)R:!0.61±0.01 K deduced in Ref.|il 
and 0, we plot Tcs vs x in Fig. We see that Tcs 
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FIG. 4: (O) Low-temperature AA(r) vs for (a) j:=0.1 and 
(b) x=0.2. The solid lines are visual aids to determining the 
range of linear fit. Tcs is defined to be the temperature where 
AA(T) starts to depart from T^-behavior. (c) (Q) Tc3{x) for 
2;=0,0.1,0.2. Line: Best linear fit to the three data points. 
Note that the line extrapolates to zero near x=0.26. 

varies linearly with x. Extrapolating the best-fit line 
yields TcS^O when a;w0.26. This implies that the low- 
T nodal phase disappears, perhaps at a quantum critical 
point, when x > 0.3, i.e. one only sees a fully-gapped be- 
havior over the whole temperature range, agreeing with 
our x>0.4 data sets. 

The continuity across the series of the first super- 
conducting transition, that we label Td, and the BCS- 
like behavior of ps over much of the T-x plane, sug- 
gest that conventional phonon-mediated superconductiv- 
ity prevails. Nonetheless, there is ample evidence for a 
second superconducting transition at Tc2 at x=0 below 
which unconventional superconductivity appears. Spe- 
cific heat measurements on Pri_j,LayOs4Sbi2 01 showed 
that the second superconducting transition at Tc2 dis- 
appears between j/=0.05 and 0.1, leaving conventional 
superconductivity for larger values of y. Figs. ISt- 



andUjD show some changes in curvature in A A close to 
T* for the x=0.1, 0.2 and 0.4 samples that could be 
indicative of Tc2, but which are not reproducible from 
sample to sample. As noted in the introductory para- 
graph, two mechanisms — spin-fluctuation and aspheri- 
cal Coulomb scattering — have been proposed to explain 
the heavy-fermion behavior and superconducting prop- 
erties of the x=0 skutterudite. One possibility is that 
the spin-fluctuation mechanism is active at high temper- 
atures where the state is thermally populated on the 
Os-rich end of the phase diagram, but is suppressed by 
decreasing temperature or as Ru doping increases the 
Fi-Ts splitting. Aspherical Coulomb scattering may re- 
main important at lower temperatures and at larger val- 
ues of X. Our data, when considered together with other 
data and theory, suggest three different superconducting 
phases: phonon-driven (conventional) across the series at 
the upper transition Td, but with spin- fluctuation and 
aspherical Coulomb scattering at the Os end giving rise 
to transitions to unconventional phases at Tc2 and Tcs- 

In conclusion, we report measurements of the 
magnetic penetration depth A in single crystals of 
Pr(Osi_2:Rux)4Sbi2 down to ~0.1 K. Both A and su- 
perfluid density ps exhibit an exponential behavior for 
the a;>0.4 samples, going from weak-coupling (a;=0.4,0.6) 
to moderate-coupling (a;=0.8). For the x<0.2 samples, 
both A and Ps vary as at low temperatures, but Ps 
is s- wave-like at intermediate to high temperatures. Our 
data are consistent with a three-phase scenario, where a 
fully-gapped phase at Td undergoes a transition to an 
unconventional phase at Tc2~Tci, then to a nodal low-T 
phase at Tcs for small values of x. The x-dependence of 
Tc3 suggests that the low-T phase disappears near x=0.3. 
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